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Abstract 



We calculate the D-wave charmonium production in e + e anni- 
hilation at BELLE and BABAR at y/s = 10.6 GeV using the non- 
' relativistic QCD factorization formalism, including color singlet and 

. color octet contributions. We analyze the contributions of various 

\ processes for 5j(J = 1,2) production. We find that both the color 

' singlet and color octet channels may give substantial contributions, 

and the production rates are estimated to be cr(5i) — 0.043 — 0.16pb 
, and 0(82) ~ 0.094 — 0.29p6, which , however, are very sensitive to 

the choice of the color octet matrix elements of the D-wave charmo- 
q h i nium states. The measurement of D-wave charmonium production at 

. BABAR and BELLE in the future will be very helpful to test the color 

<D \ octet mechanism and to determine the color octet matrix elements. 

■ PACS number(s): 12.40.Nn, 13.85.Ni, 14.40. Gx 



\ Studies of heavy quarkonium production at high energies are hot topics 

^ ' in recent years. In the conventional picture, the heavy quarkonium produc- 

tion is described in the color-singlet model|rj. In this model, it is assumed 
that the heavy quark pair must be produced in a color-singlet state at short 
distance with the same angular-momentum quantum number as the charmo- 
nium which is eventually observed. However, with the recent Tevatron data 
on high pt J production, this color-singlet picture for heavy quarkonium 
production has become questionable. The observed cross sections are larger 
than the theoretical prediction of the color-singlet model by a factor of about 
30 ~ 50@. This is called the J/ip (ip') surplus problem. On the theoretical 
side, the naive color-singlet model may be supplanted by the nonrelativis- 
tic QCD (NRQCD) factorization formalism ||, which allows the infrared 
safe calculation of inclusive charmonium production and decay rates. In this 
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approach, the production process is factorized into short and long distance 
parts, while the latter is associated with the nonperturbative matrix elements 
of four-fermion operators. So, for heavy quarkonium production, the quark- 
antiquark pair does not need to be in the color-singlet state in the short 
distance production stage, which is at the scale of 1/m.q (mg is the heavy 
quark mass). At this stage, the color configuration other than the singlet, 
i.e. the color-octet is allowed for the heavy quark pair. The later situa- 
tion for heavy quarkonium production is called the color-octet mechanism. 
In this production mechanism, heavy quark-antiquark pair is produced at 
short distances in a color-octet state, and then hadronizes into a final state 
quarkonium (physical state) nonperturbatively. With this color-octet mech- 
anism, one might explain the Tevatron data on the surplus of J/ip and if)' 
production 0. [5|. 

Aside from the J /if} production at Tevatron, various heavy quarkonium 
production processes have been studied to test the color octet mechanism. 
Among them, the charmonium production in e + e~ annihilation is particu- 
larly interesting, since in this process, the parton structure is simpler, and 
there is no higher twist effect, so the theoretical uncertainty is smaller and 
it can even be used to extract the color-octet matrix elements. J /if) pro- 
duction in e + e~ annihilation process has been investigated within the color- 
singlet model g 0, g, H, [TC|, [TTJ and color-octet model 0, |EJ 0. Recently, 
BABAR[[15] and BELLE[1(J have measured the direct J/if> production in 
e + e~ annihilation at a/s = lO.QGeV. The total cross section and the angu- 
lar distribution seem to favor the NRQCD calculation over the color-singlet 
model|l5|, but some issues still remain. The P-wave charmonium Xcj pro- 



duction in e + e~ annihilation has been discussed in |T7|]. The total x&,2 cross 
sections are dominated by color-octet process, because the C-parity sup- 
presses the process e + e~ — > ccgg, which dominates the J /if} production in 
the color singlet sector. In order to further test the color octet mechanism, 
in this paper we study the D-wave charmonium production in e + e~ annihi- 
lation. As pointed out by the authors in |L9| and p0[] , due to the color octet 
mechanism the D-wave charmonium production rates could be as large as 
J /if} or if}' at the Tevatron, and could be enhanced by two orders of magni- 
tude at the fixed target experiments. Therefore, it would also be interesting 
to study the D-wave charmonium production in e + e~ annihilation. 

We would like to mention that among the three spin-triplet D wave char- 
monium states the 3 D 2 is the most promising candidate to discover. Its mass 
is predicted to be in the range of 3.81 — 3.84GeV in the potential model calcu- 
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lation | IBfl , which is above the DD threshold but below the DD* threshold. 



The state 3 D 2 (J PC = 2 ) is forbidden to decay into DD because of the 
parity conservation. So this state is predicted to have a narrow decay width 
of 300 — AOOkeV, and we may easily tag this state through the dominant 
decay channels, e.g. the El transitions into XcJ=i 2 states, and the hadronic 



transition into J /if) + 7T7r(see [19] and [2C] for detailed discussions). For the 
3 Di state ^>"(3770), it is just above the DD threshold and thus has a quite 
narrow width of about 23MeV and may be tagged via if>" ~^ DD. As for 
the 3 D 3 (J PC = 3—) state, its mass is far above the allowed decay channel 
DD threshold and it has more decay modes, therefore it is expected to be a 
wide resonance and difficult to observe. So, phenomenologically we are only 
interested in the 3 D 2 as well as 3 Di states, which will be discussed below. 

The power of the NRQCD formalism stems from the fact that factoriza- 
tion formulae for observables are expansions in the small parameter v ,where 
v is the average relative velocity of the heavy quark and antiquark in quarko- 
nium bound state. For charmonium v 2 ~ 0.3 , and for bottomonium v 2 ~ 0.1 



. NRQCD velocity-scaling rules allow us to estimate the relative size of 
various NRQCD matrix elements. This information, along with the depen- 
dence of the short-distance coefficients on a s , permits us to decide which 
terms must be retained in expressions for observables to reach a given level 
of accuracy. At low orders, factorization formulas involve only a few matrix 
elements, so several observables can be related by a small set of parameters. 

In NRQCD the Fock state expansion for the physical D-wave charmonium 
(denoted by Sj) is 

\Sj) =0(l)\cc( 3 Dj,l)) 
+0{v)\ce{ 3 Pj,,8)g) 

+0(v 2 )\cc( 3 S 1 ,8orl)gg) + ---. (1) 

The striking feature of this expansion is that although the amplitudes of 
these Fock states in the expansion are different in v the contributions to the 
production rates of all these terms are essentially of the same order of v 7 . 
This is very different from the case of J /if) and if)' '. 

According to the NRQCD factorization formalism ||], the production pro- 
cess e + e~ — > 5 j + X can be expressed as the following form, 

da(e + e- ^ 5j + X) =J2F(e + e- ^ n + X)(O s n J ). (2) 
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Here, n denotes the cc pair configuration in the expansion terms of Eq. ([I]) 
(including angular momentum 2S+1 Lj and color index 1 or 8). F(e + e~ — > n + 
X) is the short distance coefficient for the subprocess e + e~ — > n + X. (O s n T ) 
is the long distance non-perturbative matrix element which represents the 
probability of the cc pair in n configuration evolving into the physical state 8j. 
The short distance coefficient F can be calculated by using perturbative QCD 
in powers of a s . The long distance matrix elements are still not available 
from the first principles at present. However, the relative importance of the 
contributions from different terms in Eq. (Q) can be estimated by using the 
NRQCD velocity scaling rules. 

The main Feynman diagrams for the production of 8j in e + e~ annihilation 
are shown in Fig.l, where (a) and (b) are the color singlet processes, and (c) 
is the color octet process. 

For the color-singlet process in Fig. 1(a), 

e + e- ^7* ^cc[ 3 Dj,l}+g + g, (3) 

With the spin-triplet case where J = 1, 2, 3, we use explicit Clebsch- 
Gordan coefficients, and get the following relations for the three cases. 

V/H lUUXW - r 3 U/2r/„ PaPp s {J z ) , / PpPp \ <J Z ) 

}^{U z \ls2m)e af3 ey - -[—J [(g ap - jfipfy + » ~ J^p) 6 " 

STY}/ C C 

g<2J,|lrfro>e$e« = ^tf^r' + ^W^V*'), (5) 
J2(SJ z \ls2m)e^e^ = e<#. (6) 

sm 

Here, e a , e a p, e a p p are J = 1, 2, 3 polarization tensors which obey the projec- 
tion relations 

EJm) (m) _ ( , PaPP \ _ -p / ? \ 

m ^"■c 

5Z e i™ )e i™/3' = ^[Paa'Vpp' +V a p'V/3a'} - l{P a pVa'l3', (8) 
rn Z 6 

^ e jrjp e L^9V = -^Pa.al'PppVppi +V aa 'Vf3 p 'V l 3 p ' + V a /3'Vf3 a 'V pp ' 



+ Vafi'Vpp'Vpa' + Vap'Vpp'Vpa' + V a p'V p a 'V p p') 
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Figure 1: The main Feynman diagrams for the production of 5j in e + e 
annihilation 



5 



— -^{VapPpaiVpipi + V a pV p p'V a ' p ' + V a pV pp 'V a ' p> 

+ V ap Vp a 'Vp' p ' + VapVpp'Ta'pi + V ap Vpp'V a 'p' 

+ VppVoca.iV pi p > + Vp p V a p'V a 'p' + VppV a p'V a 'p'). (9) 

We choose the parameters: 

m c = 1.5GeV,a s (2m c ) = 0.26, 

\R'^(0)\ 2 = 0M5GeV 7 (10) 

For the color-singlet matrix element, we use 

(Oi J CDj)) = 15(2J 4 + 1)iYc |^(0)| 2 . (11) 

Because these color-singlet processes have the infrared divergence in- 
volved, we introduce an infrared cutoff, which can be set to m c on the energy 
of the outgoing gluons in the quarkonium rest frame. For these color-singlet 
contributions, we finally obtain 

<?singlet{5l) = 0.027pb 

Vsinglet(5 2 ) = 0M7 P b (12) 

For the color singlet process in Fig. 1(b), we have estimated it in the frag- 
mentation limit, and found that its contribution to the S\ production cross 
section is smaller by a factor of about 60 than that for the J/ip production 
cross section which is about 0.1 ~ 0.2pb. This is consistent with the result 
given in|22|]. So the contribution of Fig. 1(b) to the 5\ production cross section 
is only 0.002 ~ 0.004p6, to 82 is 0.003 ~ 0.006p6, therefore is negligible. 

We now discuss the color octet processes in Fig. 1(c) 

e+e- ^9 + cc{ 2S+1 Lj,8], (13) 

we readily have 

a ( e + e - ^5j + g) = C s (0 5 8 > CS )) + C P (O 5 8 '( 3 P )), (14) 

with 

c = wj^izr, (15) 

3 s z m 



2567r 2 e c 2 « 2 a sr (l-3r) 2 6(1 + 3) 2(1 + 3r + 6r 2 , 

C p = 7TY~3 Hj + + 1 h ( 16 ) 

9s^m d 1 — r 1 — r 1 
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where r = m 2 /s, m is the mass of 5j, and s is the e + e cm. energy squared. 
Here we have used the approximate heavy quark spin symmetry relations 

{O 5 8 '( 3 Pj))^(2J+l)(O 5 8 '( 3 P )) (17) 

Then we calculated the cross sections of Sj production. The cross section is 

a(e + e' ^5j + X) = Co(0%CS 1 )) + C 1 (0 5 8 '( 3 P )) + C 2 (O s 1 '( 3 Dj)) (18) 

Using the NRQCD velocity scaling rules, we roughly estimate the 5j color 
octet matrix elements by relating them to the ip' color octet matrix elements: 

(O s 8 J ( 3 S 1 )) « (2J _ + 1 W t 3 ^)) ~ (2J _ + 1} x 4.6 x 10- 3 GeV 3 (19) 
o o 

(O^( 3 P )> » (2J r +1) (O 8 ^( 3 P )) ^ (2J r + 1} x 1.3 x 10- 3 Ge\/ 5 (20) 
5 5 

Using those color-octet matrix elements, we get 

(T^OJi) ~ 0.016p6 (21) 
<w(5 2 ) ^ 0.027pb (22) 

Then we get the total cross section: 

a(e + e~ -> 7* -> tfj. + X) ~ 0.043pfe 

(j(e + e- -> 7* -> 5 2 + X) ~ 0.094p6 (23) 

However, if we use a more radical choice of the 5j color octet matrix 
elements by relating them to the J / ip color octet matrix elements, we would 
get much larger values for the 5j production cross sections, cr(5\) — 0.16p6 
and a(S 2 ) ^ 0.29pb, because the J ftp color octet matrix elements are much 
larger than the ip' color octet matrix element s((O 8 ^( 3 P )} ~ 1.1 x lCT 2 GeU 5 , 
elements((Og ( 3 -Po)) — 1-3 x 10~ 3 GeU 5 ). Therefore the measurement of Sj 
production rates in the future will be very helpful to determine the values of 
the Sj color octet matrix elements. 

In summary, we have estimated the D-wave charmonium 5j(J = 1,2) 
production rates in e + e" annihilation at y/s = 10.6 GeV. We find that both 
color singlet and color octet may give substantial contributions to the pro- 
duction cross sections, of which the values are sensitive to the choice of the 
5 j color octet matrix elements. We hope that these results can be tested 
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with higher statistic e + e~ annihilation data in the future at BABAR and 
BELLE. For 5 2 , its branching fraction of decay mode J/-ip7T + 7T~ is estimated 
to be B(82 — > J/^7T + 7r _ ) ~ 0.12^3], which is only smaller than that of 
B{ip' — > J/ipir + ir~) = 0.324 ± 0.026 by a factor of 3. We compare the pre- 
dicted production rate of 2 D-wave charmonium with that of ip'. With 
the integrated luminosities of about 30/6 -1 at ^/s = 10.6 GeV, BELLE 
gives a(e + e~ — > ip' + X) ~ 0.67 ± 0.09p6 with 143 ± 19 ip' events decay- 
ing to J/ip7i + 'K~ fllEf . As a rough estimate, if we choose (7(82) = 0.094p6, 
about 7 events of 2 state decaying to J/ipn + n^ will be detected, and for 
(7(82) = 0.29pb (with larger values for the color-octet matrix elements), it is 
about 23 events. With more data available at the B Factories in the future, 
it will be possible to detect this 2 D-wave charmonium state (especially 
in the latter case). The 3 Di cc state ?/> // (3770) could also be detected via 
ip" -> DD decay. 
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